AD  744674 


NIL  Report  5259 


ON  THE  THEORY  OF  FLEXURAL  PiEZOCERAMIC 
CIRCULAR  PLATE  SOUND  RADIATORS 


S.  Hanish 

Transducer  Branch 
Sound  Division 


February  26,  1959 


C.  S.  NAVAL  RESEARCH  LABORATORY 
WasAiaptOft,  D.C. 

N’.VT  TEChNiO A!. 

!Nt-  v  •’»?;  ON  SE^X  iCE 


approves  ?c?  r.T" 
dzsirif-uiIvS  •  ■- 


CONTENTS 


Abstract  ii 

Proolem  Stab-.s  ii 

Authorization  ii 

INTRODUCTION  1 

CENTRALLY  SUPPORTED  FLEXURAL  DISKS  1 

The  Structure  1 

Equations  of  State  2 

Coefficient  of  Electromechanical  Coupling  4 

Mechanical  Displacement  Under  Forced  Electrical  Drive  4 

Transduction  Ratio  9 

Radiated  Acoustic  Power,  Mechanical  Reactive  Power,  Kinetic  Energy, 

Resonant  Frequency,  Acoustic  Pr*  ;ure  and  Mechanical  q  of  a 
Centrally  Supported  Disk  Whose  Deflection  Curve  is  Parabolic 
(Second  Order  in  r).  Infinite  Baffie  Present.  12 

Mechanical  Displacement  Under  Forced  Mechanical  Drive  and  Receiving 

Response  15 

Formulas  Describing  the  Acoustic  Performance  of  a  Metal-Ceramic 
Bilamellate  at  Mechanical  Resonance,  '  dialing  Sound  into  a  Liquid 
Medium.  Parabolic  Deflection  Curve  Assumed  and  Infinite  Baffle 
Present.  16 

Conclusion  17 

EDGE-SUPPORTED  BILAMELLATE  DISKS  18 

Coefficient  of  Electromechanical  Coupling,  Bending  Moments,  Dynamic 

Equation,  ana  Solutions  18 

Acoustic  Power,  Mechanical  Reactive  Power,  Kinetic  Energy,  Mechanical 
Resonant  Frequency  and  Mechanical  q  of  a  Simply  Supported  Disk  Whose 
Deflection  Curve  is  Parabolic  (Second  Order  in  r).  Infinite  Stiff  Baffle 
Present.  20 

Formulas  Describing  the  Acoustic  Performance  of  a  Simply  Supported 
Ceramic  Bilamellate,  at  Mechanical  Resonance,  Radiating  Sound  into  a 
Liquid  Medium.  Parabolic  Deflection  Curve  Assumed,  and  Infinite 
Baffle  Present.  22 

Conclusion  23 

LIST  OF  SYMBOLS  24 

REFERENCES  26 


i 


ON  THE  THEORY  OF  FLEXURAL  PIEZOCERAMIC 
CIRCULAR  PLATE  SOUND  RADIATORS 


INTRODUCTION 

A  convenient  source  of  underwater  sound  in  the  range  of  100  cps  to  50  kc  iu  the 
family  of  flexural,  circular,  piezoceramic  bilamellate  disks.  These,  in  diameters  of  0.5 
to  25  cm  can  radiate  up  to  1/2  to  1  watt  per  square  centimeter  for  applied  electric  fields 
of  150  to  250  volts  rms  per  millimeter  of  half  thickness,  pulsed  power.  A  simplified 
theory  of  the  performance  characteristics  of  such  sources  is  presented  in  this  report. 

In  the  first  part  of  the  theory  the  piezoceramic  disk  is  centrally  supported,  and  in  the 
second  part  the  disk  is  simply  supported  on  its  edge.  Both  cases  are  treated  only  with 
reference  to  the  presence  of  an  Infinite  stiff  baffle  enclosing  the  half  space  into  which  the 
sources  are  radiating  acoustic  energy.  The  analysis  proceeds  on  the  assumptions  of  a 
Class  I  (Fischer's  classification)  piezoelectric  transducer  and  not  on  the  basic  laws  of 
the  Class  HI  transducer,  to  which  the  electro strictive  types  belong.  This  has  been  done 
to  avoid  transduction  ratios  and  lumped  parameters  which  are  frequency  dependent.  As 
a  result  of  this  choice,  the  material  constants  appearing  in  the  course  of  the  paper  are 
only  "effective"  piezoelectric  moduli,  to  be  measured  on  particular  samples  of  polarized 
ceramics  by  techniques  conventionally  applied  to  true  piezoelectric  crystals.  In  short, 
the  ceramics  have  been  converted  to  "effective"  crystals,  and  standard  piezoelectric 
theory  has  been  applied  to  their  behavior. 


CENTRALLY  SUPPORTED  FLEXURAL  DISKS 
The  Structure 

The  transducer  (Fig.  1)  is  a  bilamellar  disk,  one  or  both  halves  of  which  may  be 
actively  polarized  electro  strictive  ceramics.  While  it  is  always  advantageous  to  have  a 
disk  of  two  active  halves,  the  use  of  a  central  support  may  require  half  of  the  bilamellar 
structure  to  be  a  passive  metal,  e.g.,  brass.  Each  face  of  an  active  ceramic  plate  is 
completely  electroded  with  fired  silver  paste.  Permanent  polarization  of  the  ceramic  iB 
in  the  z  direction,  normal  to  the  plate.  The  structure  is  driven  into  forced  flexural  vibra¬ 
tion  by  the  application  of  an  alternating  electric  field  E  applied  across  the  thickness  of 
the  ceramic  plate.*  Upon  submersion,  the  disk  radiates  sound  into  a  semi-infinite  liquid 
medium  of  characteristic  impedance  pm cw,  from  a  circular  hole  of  equal  diameter  in  an 
infinite  stiff  baffle. 

Figure  2  shows  the  conventions  of  signs  which  hold  in  this  analyst.  The  neutral 
axis  of  bending  is  set  at  z  =  o.  For  small  displacements  w  in  the  positive  z  direction,  the 
displacement  u  in  the  x  direction  is  u=  -zOw/Bx),  where  3w/Bx  is  positive,  as  shown.  The 
strain  Sxx  in  the  x  direction  is  Sxx  =  3u/3x  =  -z02w/bx2) .  We  note  that  b2w/bx2  is  nega¬ 
tive  when  w  is  positive  and  that  the  strain  is  positive  (tension)  for  positive  z.  Similarly, 
for  displacement  v  in  the  y  direction,  the  strain  in  the  y  direction  is  Syy  =  -zOv/b y)  = 
-z(32w/3y2).  A  positive  bending  moment  Mx  or  My  is  one  which  causes  a  deflection  of 
the  disk  in  the  positive  z  direction,  i.e.,  downward,  as  shown,  and  induces  resisting 
stresses  Txx  or  Tyy . 


♦See  List  of  Symbols  at  end  of  report. 
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Figure  1  -  Schematic  cross  section  of  a  piezoceramic -brass 
circular  plate  sound  radiator 


Figure  2  -  Diagram  illustrating  deflection 
parameters  cr  an  elemental  section  of  the 
circular  plate  sound  radiator 


The  same  conventions  hold  for  a 
cylindrical  coordinate  system  in  which 
the  origin  corresponds  to  the  origin  of 
the  xyx  Cartesian  system. 


Equations  of  State 

The  circular  bilamellar  disk  (half 
thickness  h)  is  assumed  "very  thin,  ''  td 
the  stress  system  correspondingly  two- 
dimensional.  This  restriction  to  stresses 
in  a  plane  is  indicated  at  any  stage  in  the 
analysis  by  a  bar  over  all  material  con¬ 
stants  (i.e.,  crr,  ztf,  etc.).  A  cylindri¬ 
cal  coordinate  system  is  used,  and,  for 
convenience  in  writing,  r,  *  are  desig¬ 
nated  by  the  numbers  1,  2,  and  3.  In  such 
a  planar  cylindrical  system  the  equations 
of  state  of  the  bilamellar  disk  under  sym¬ 
metrical  electrical  excitation  may  be 
immediately  taken  from  the  tensor  form 
used  by  Mason  (1). 
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The  negative  sign  appears  in  the  term  -e31  C3  for  the  following  reason.  H  a  cylindri¬ 
cal  pill  of  material  is  scooped  from  the  interior  of  the  ceramic,  and  if  a  positive  gradient 
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e3  is  applied  across  the  thickness,  the  pill  is  assumed  to  expand  in  the  *  direction,  the 
polarization  being  correctly  oriented  for  this  expansion.  Now  if  the  material  adjacent  to 
the  pill  in  the  z  direction  is  hindered  from  moving  (condition  of  zero  strain),  it  will 
develop  a  resisting  compressive  (negative)  stress.  Similarly,  simultaneous  contraction 
of  the  pill  in  the  r  and  5  directions  will  induce  tension  stresses  in  the  adjacent  material 
which  are  the  resisting  stresses  Tj  and  T2.  Since  the  "effective"  piezo  modulus  e31  is  a 
negative  number,  a  minus  sign  is  added  to  the  term  e31  E3  so  as  to  render  Tt  and  T2 
positive. 


Corresponding  to  the  internal  resisting  stresses  and  T2  due  to  the  applied  field  E3 
are  the  internal  bending  moments  Mr,  V*  resulting  from  bilamellar  action.  By  suitable 
integrations  over  the  thickness  of  the  plates,  expressions  are  obtained  for  these  bending 
moments  in  terms  of  displacement  derivatives. 


z  At 


T2  i  dr 
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(*> 
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rf»\  _  g3i 
dr  j  2 


1  d»\  *31h~E3' 
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where 

Mr  is  the  internal  sending  moment  per  unit  of  circumferential  length 

M*  is  the  internal  bending  moment  per  unit  of  radial  length 

D  is  the  flexui  al  constant,  -  ( 2  '31  h3  c  E, 

is  Poisson's  ratio,  which  for  planar  stress  is  c,E2  c/j 

t  is  1  or  2,  dc.*  .-Ang,  respectively,  on  whether  the  backing  plate  is  inactive  (i.e., 
is  metal)  or  a  r.  -  (i.e.,  is  a  ceramic  plate  polarized  for  bilamellar  action), 
assuming,  if  ao/.-.  that  active  halves  are  connected  in  paralleL 


Similarly,  the  average  A-  trie  displacement  in  one  plate  is 


e3»  h  /  d2w 
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(3) 


As  the  radiation  ■■■■:  induces  a  real  surface  pressure  on  the  face  of  ti.e  disk 

exposed  to  the  medium,  an  expression  is  needed  lor  the  internal  resisting  shear  or  per 
unit  of  circumferential  length.  Since  the  applied  electric  field  E3  is  independent  of  the 
radial  coordinate  r ,  we  may  use  conventional  thin-plate  theory  and  write  (2) 
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dr* 
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(4) 


Equations  1  to  4  are  the  foundation  equations  upon  which  the  frame  of  the  analysis  is 
built.  They  carry  with  them  the  limitations  that  restrict  the  final  formulas  to  narrow 
grounds  of  validity.  In  particular,  they  are  piezoelectric  in  origin,  linear,  and  planar, 
the  electric  fields  are  low  in  frequency  and  small  in  magnitude,  the  sound  fields  are  due 
to  infinitesimal  displacements,  an  infinite  stiff  baffle  is  present,  and  edge  effects,  cement 
between  plates,  etc.,  are  considered  to  affect  the  results  in  minor  ways  only. 
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Coefficient  of  Electromechanical  Coupling 

The  mechanical  and  electrical  equations  (Eqs.  2  and  3)  for  the  state  of  flexure  must 
be  satisfied  together;  that  is,  they  arise  from  independent  laws  of  nature  and  are  not 
reducible  to  a  single  statement  except  by  simultaneous  solution.  If  the  expressions  for 
Mr  and  are  added  and  the  resultant  solved  for  the  displacement  derivatives,  we  obtain 


d^r  1  dw  .  Mr  *  Mf  *31 E3  i 
r  dr  '  Dfl**)  D(l+v)  ’ 


i.e.,  mechanical  curvature  of  the  disk  is  expressed  in  terms  of  the  internal  bending 
moments  representing  the  mechanical  and  electrical  fields.  Since  the  same  quantity  on 
the  left  appears  in  the  electi  ic  field  equation  (Eq.  3),  we  have,  upon  substitution,  a  single 
equation  combining  the  influence  of  both  fields. 


D3 


2  \D 


+  ( 


»*3 
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fi 
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(5) 


This  equation  is  a  statement  that  a  charge  per  unit  area  on  the  electrodes  of  one 
plate  accr  mp antes  the  condition  of  flexure.  The  first  term  on  the  right-hand  side  is  the 
mecbani.al  or  motional  charge,  and  the  second  term  is  the  charge  due  to  applied  field. 
The  coup’5""  effect  of  the  mechanical  field  appears  as  an  increment  in  the  clamped 
dielectric  cons.-oit  *33»  which  takes  the  form  of  the  second  factor  in  the  enclosed  term. 
Such  a  group  of  symbols  appears  in  planar  piezoelectric  analysis  and  i«  cr:.»emently 
designated  ar  th'  coefficient  of  electromechanical  coupling  in  the  flexure  mode.  Bech- 
mann  (3)  defines  the  mixed  planar  coefficient  of  electromechanical  coupling  (vp2  )miM  by 
the  relation 


<W2> 


2e 


-2 


31 


P  mi>  *U<l  +  v'>tU 

For  the  state  of  flexure,  therefore,  the  coefficient  of  electromechanical  coupling  k{2  is 

k'2  (6) 


We  note  that  this  equation  is  independent  of  the  manner  in  which  the  disk  is  mechani¬ 
cally  supported.  When  the  structure  has  one  active  half  and  is  metal  backed,  kf2  = 

3/8  (kJ')mix  .  Equation  (6)  reveals  an  important  feature  of  flexural  sound  radiators, 
namely  that  the  interconversion  of  electrical  and  mechanical  energy  is  markedly  reduced 
by  flexure  from  the  latent  potential  of  the  planar  state  to  3/8  ;  of  this  potential.  This 
reduction  of  coupling  is  balanced,  of  course,  by  the  increased  compliance  of  mechanical 

structure  and  the  consequent  lowering  of  natural 
frequencies  of  vibration.  In  many  applications, 
the  latter  advantage  more  than  overrides  the 
resulting  poorer  power-handling  capacity  of  the 
structure. 

Mechanical  Displacement  Under 
Forced  Electrical  Drive 

The  mf  chanical  displacements  of  a  piezo- 
ceramic  disk  which  radiates  sound  into  a  semi¬ 
infinite  liquid  medium  may  be  found  by  applica¬ 
tion  of  conventional  elastic  and  acoustic  theory. 
Figure  3  shows  an  elemental  volume  of  the  disk 


z 


Figure  3  -  Elements',  volume 
of  a  circular  plcte  sound 
radiator 
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in  cylindrical  coordinates;  m  is  the  mass  per  unit  of  area  and  the  total  differential 
bending  moment  per  unit  of  circumferential  length.  Equating  internal  force  per  unit  cir- 
cumferential  length  on  n  (ft  in  the  x  direction)  to  the  acceleration  of  n ,  we  have 

dSL. 

F T  =  -  =  nfrd^dr,  2-|. 

1  df  dt2 

Now  to  quantities  of  the  first  order  in  infinitesimals,  we  see  from  Fig.  3  that 

«,*[(«,  ♦  ^  dr)  (r  ♦  dr)  AS  -  H/tAff)  -  2  ^dr  ^')+  Qr  rdrdsj 

(  HM  \ 

=  +  r  -  H#  +  rOr  j  Hr  AS. 

Letting  the  primed  symbol  represent  differentiation  with  respect  to  r,  the  force- 
acceleration  equation  becomes 


2m;  tm;  *,q;  +  Qr  =  _ 

Although  m'  and  m‘  are  independent  of  the  applied  electric  field,  the  magnitude  of  the 
elastic  constants  cn ,  c12,  etc.,  depends  upon  this  field.  For  low  electric  fields,  we 
assume  D3  very  small  and  use  the  constants  eft ,  dft  instead  of  the  constant  eft,  z12, 
etc.,  win  re  by  definition 

cf,  =  ?ft  (1  *  k,2) 
ei°i  =  ei2  «  tkf2>- 

For  the  state  of  larger  electric  fields,  03  is  hardly  zero,  and  we  must  use  a  value  of 
cjj  that  lies  between  eft  and  eft.  In  all  cases,  however,  k2  is,  at  most,  about  0.10,  and 
the  correct  choice  for  (he  constants  is  seen  to  be  no  critical  matter.  Substituting  for  il'. 
Up,  thus  modified  for  the  particular  magnitude  of  drive,  we  obtain 

d4w  ^2  d3«  _  1  d2w  +  J_  dw  +  _m  d2*.  Qr  „  1  d<?r 

dr4  r  dr3  r2  dr2  r3  ^  D*  dt.2  rD*  D* 


where 


_  _  2.  3?d 
—  n  c_ .  « 
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The  internal  shear  force  Qr  per  unit  of  circumferential  length  depends  upon  q(r),  the 
external  load  intensity  per  unit  area.  Static  ec(uilibrium  requires  that 


2rrQr 


2nrAr  q(r)  . 


from  wh'~h 


3l  +±  dQr  -  q(r) 

rD*  D*  ’  D* 

From  the  assumptions  of  thin-plale  theory,  we  restrict  q(r)  to  acoustic  loads;  i.e., 
qfr)  is  the  reaction  acoustic  pressure  p(r>.  In  the  condition  of  steady-state  sinusoidal 
vibration,  the  equation  of  motion  becomes 
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'(a2_ 

\dr2 
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d_ 

dr 


+  k  2 M  tf  -  P^f  ,c“l 

vj  F  ’ 


(7) 


where 


3pp  ‘-■2 

hT^lT- 
n  Ji 


The  surface  pressure  p(t,o>)  may  be  determined  from  the  following  considerations. 
Since  each  element  of  area  d  A  and  veloc**/  w  radiates  a  spherical  wave  into  the  half 
space,  the  surface  acoustic  prepare  is  a  summation  of  the  pressure  interaction  effects 
of  these  waves.  We  have,  from  *  is  rher  (4), 


p(r,t) 


2rr 


w 


dA, 


where  R  is  the  distance  from  an  arb:  cary  point  P  in  the  plane  of  the  iisk  to  an  element  of 
area  RdRd<p  and  <?  is  the  angle  between  R  and  the  line  P  drawn  to  the  origin  of  coordinates. 
For  the  steady  state,  if  is  equal  to  j&*W,  and  sinct  p(r,u)  is  a  negative  quantity,  we  obtain 
the  final  e*,  stic-acoustic  equation  motion 


,4_  '  _  -ikS 

*  •  kj  V  :  — - — —z—  |  W  e  d  H  dd> , 


2tr  D* 


where  v"w  -  !-4  w  equate.  -hand  side  of  Ec  7. 


(8) 


Except  for  very  sin.*  :e  geometries  and  for  simple  displacement  distributions,  this 
integro-differ  -Jtial  equation  is  intractable.  We  can  avoid  this  difficult  mathematical 
situation  by  making  the  right-hand  side  independent  oi  the  variable  r,  which  means  in 
effect  independe-.  -f  the  variable  W.  The  reaction  pressure,  which  is  known  to  be  a  func¬ 
tion  of  the  coordinate  r  and  the  frequency  u,  is  replaced  by  an  approximate  pressure 
constant  with  radius  but  strongly  dependent  on  the  wavelength  of  radiation,  that  is,  on 
frequency.  Instead  then  ol  seeking  an  exact  solution  to  the  radiation  problem  posed  by 
Eq.  (8),  we  will  seek  the  solution  to  a  germane  problem  in  which  the  external  reaction 
pressure  of  the  liquid  medium  is  arbitrarily  defined  to  be  independent  of  r,  while  still 
retaining  dependence  on  «.  In  place  of  Eq.  (7),  therefore,  we  will  set  the  equation  of 
motion  to  be 


V4f  - 


k  *TI  = 


p(a.oQ 

D* 


where  a  is  the  outside  radius  of  the  disk. 


(9) 


The  reaction  pressure  pfa.u)  can  be  only  an  average  estimate  of  the  true  pressure 
p(r.<o).'  A  convenient  procedure  for  determining  p(a,&>)  is  to  assume  a  displacement 
curve,  V(r)  for  the  flexed  disk,  conformable  to  the  boundary  conditions,  and  perform  the 
integration  required  by  the  right-hand  side  of  Eq.  (8).  The  result  of  this  first  step  is  the 
pressure  function  p(r.w).  _  Me  Lachlan  (5)  carried  through  such  an  integration  for  the 
assumed  distribution  w  =  wo[i  -  r(r2/a2)]  and  found  that 


=  p-c-w 


Ld 


,J  2S 

+  j  1**1  "  3f  *3  J  57  «s 


-)]• 


(10) 


where  g2  and  g6  are  known  algebraic  functions  of  the  radius  r,  and  cx  and  gs  are  known 
hypergeometric  functions  F(a,£,y,6)  of  the  same  coordinate.  The  symbol  *  stands  for  the 
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wave  parameter  k».  The  infinite  series  of  Eq.  (10)  is  applicable  to  an y  size  disk  having 
the  velocity  distribution  prescribed.  In  general,  however,  the  frequencies  corresponding 
to  the  grave  mode  of  vibration  are  so  low  that  k«  is  less  than  1/2.  We  assume,  then,  that 
ka  «  1/2  and  write,  after  Me  Lachlan  (5), 


P(r.t)  =  P„cwW0  +  j(ka)gtJ  (11) 

for  the  case  of  the  centrally  supported  disk.  The  exact  expression  for  gt  is 


The  average  pressure  p(a,t)  is  found  by  averaging  Eq.  (11)  with  respect  to  disk 
area.  The  result  is 


P'a.t)  =  p.ewt0  ♦  j  ^  k»]  • 


(12) 


A  convenient  symbol  for  the  impedance  portion  of  this  equation  is  Z(a.w),  defined  by 
the  relation 

p(a.t)  =  Z(*,oS)  * 


Z(a.*>)  =  pvcw  +  j  ^  kaj  . 


Having  now  defined  the  nature  of  the  expression  p/a.t),  we  return  to  Eq.  (9)  and  sub¬ 
stitute  our  results,  recalling  that  in  the  steady  state,  V  =  ju«o.  We  have,  then. 


V>  -  k<W  = 

D* 


(13) 


This  is  an  inhomogeneous  differential  equation,  the  right-hand  side  of  which  is  independ¬ 
ent  of  the  coordinate  r .  A  particular  solution  is  easily  seen  to  be 

j&>W0Z(a,a») 

D’k«  ' 


The  complementary  solution  of  the  homogeneous  equation  is  found  in  various  texts, 
e.g..  Ref.  6.  Adding  the  particular  solution  to  the  complementary  solution,  we  obtain  the 
complete  solution 


W(r,*)  =  aj0(kr)  +  0Yo(kr)  ♦  yl0(kr)  *  SK^kr)  - 

where  a ,  /3 ,  y ,  and  i  are  four  constants,  to  be  determined  from  the  conditions  of  support. 

The  boundary  conditions  for  the  case  of  a  flexural  disk  supported  by  a  rigid  built-in 
central  pin  (radius  b)  are  easy  to  formulate,  but  they  may  not  correspond  to  the  actual 
stress  configuration  because  of  local  yielding,  rotary  inertia  effects,  etc.  At  the  risk  then 
of  stipulating  what  may  be  an  ideal  situation  at  the  point  of  support,  we  specify  that 
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at  r  =  b  at  r  =  a 

1.  W  =  0  3.  M_  =  0 

2.  =  o  4.  Qr  =  0  . 

Four  conditions  are  thus  available  for  determining  the  four  constants  a,  s,  y,  and  S 
by  employing  explicitly  the  general  solution  contained  in  Eq.  (14).  Upon  performing  the 
necessary  derivations,  we  come  to 

jaiZ(  a,u>)  W 

aj0(k,b)  +  /SYoOcjb)  +  ylofkjb)  +  6K0(kjb)  =  - j  —  (15a) 

“k(ajj(k jb)  +  iYjfkjb)  -  yl^kjb)  -  iK,(k,b)]  =  0  (15b) 


aAj  +  *k2  +  rA3  6A4  =  — f- 

*11  kl  h 


where 


k  3  [aBj  +  3B2  +  yB3  ♦  6B4]  =  0 


(I5d) 


Aj  =  -Jo(kia)  *  (1  -  vlJjfk^yki* 

A2  =  -Y0(kj*)  +  (1  -  vlYjCkj.lAi* 

A3  =  -IoO^a)  -  (1  -  v)X1(kt«Vk1« 

A4  =  K^kja)  +  (1  -  s/lKjfkjalAja 

Bj  —  J|(kj«),  B2  -  Yj(kjs),  B3  -  I j(tja),  B4  —  -Kj('k ^8;  ■ 

Simultaneous  solution  of  Eqs.  (15a,  b,  c,  d)  is  rapidly  performed  by  use  of  Cramer's 
rule.  We  can,  however,  simplify  the  results  by  noting  that  for  centrally  supported  disks 
the  ratio  b/a  is  considerably  less  than  unity,  so  much  less,  in  point  of  construction,  that 
we  may  assume  that  b— *0  without  major  error.  As  a  consequence  of  this  choice,  Eq. 
(15b)  reduces  tc 


or 


LoYj(k,b)  +  SKjtkjb))  =  0. 


Substituting  tnis  result  in  Eq.  (15a),  and  again  letting  b— *o,  we  obtain 

j  wZ( 

a  +  v  - - . 

D*k« 

Completing  the  solution  for  the  two  remaining  conditions,  Eqs.  (15c)  and  (15d),  and 
substituting  the  results  in  Eq.  (14),  we  obtain 
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~n  (B2  ♦  I  B4)  a  [b3  (*a  +  |a4)  -  a3  (b2 


[J0(kir)  -  I0(k,r)]- 


A3(B,  -  B3)  -  B3(A,  - 


2  1 
Y0(k,r)  +  ±  K0(k1r )j 


where 


+  0[Io(k1r)  -  1]-, 


-  2  e31E3 

V  “  «  cD  hk  2 
11  1 


jo>Z(*,o>)f 


A  =  (A,  -  A3)  (b,  *|b4)  -  (Bj  -  B3)  (aj  ♦!*)• 

The  important  quantity,  of  course,  is  the  displacement  at  the  outer  edge,  wo.  At .  =  a, 

V  _  ~v  {(B2  B«)  ~  *o(ki«>I  ~  <Bi  ~  1;3>  [vki«)  +  'f  Ko(kia)]} 

A  |B3  [a2  + 1  a4]  -  A3  [bj  ♦ |  b4]  tj0(k,.)  -  i0(klB)] 

-  {A3<B1  -  B3)  -  B3<Ai  -  A3)J  [YoCkja)  ♦|K0(kta)] 

A  j  «Z(a,w)  rr  /t.  -x  ,i 

- 55-j-ji -  ll0(kl«)  ~  H  •  (17) 

The  edge  (maximum)  displacement  wo  is  proportional  to  the  electric  gradient  e3  (that 
is,  proportional  to  the  constant  v).  In  the  absence  of  acoustic  load  (i.e.,  when  Z  =  0)  the 
expression  for  A  becomes  zero  for  an  infinite  number  of  values  of  ka.  The  lowest  of 
these  (excluding  ka  -  0)  corresponds  to  the  grave  (or  umbrella)  mode  of  vibration;  i-.e 
remaining  values  of  ka  correspond  to  modes  of  vibration  consisting  of  a  successively 
increasing  number  of  nodal  circles.  With  the  accession  of  an  acoustic  load  the  denomi¬ 
nator,  for  certain  values  of  ka,  reaches  a  minimum,  but  it  may  never  be  zero,  since 
Z(a,&>)  is  a  complex  quantity.  We  note  in  particular  that  it  is  the  resistive  part  of  Z(a,&>) 
which  contributes  an  imaginary  term  to  the  denominator  and  diat  a  purely  reactive  load 
will  not  restrain  the  motion  of  the  disk  at  mechanical  resonance. 


Transduction  Ratio 

The  time-varying  charge  Q  accumulating  on  the  electrodes  of  one  plate  is  obtained 
by  integrating  the  expression  for  the  dielectric  displacement,  Eq.  (3),  with  respect  to 
plate  area.  We  have  (assuming  b— *0) 


Q=j[  D3  dA  =  -e31  hn  |  +  7  3?)  rdl 

Q  =  -tre31h.  [gj^  ♦7T«2e$3E3. 


+  TT12  €„E, 


+  7T*^  CJ3  E3  • 


Trr 


"w 
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The  alternating  current  due  to  an  applied  voltage  V3  (across  one  plate)  is  therefore 
f-i?  (b2  +  |  B4) 

i  =  j  o.77e31hk1a  •( - » - r - (Jjfkja)  +  Ij(kji 


,a» 


l(Bl  -  B3) 


[Yjfkja)  +  f  K,(ka)] 


*  [°3  ^  *  » a<^  *»  k 4  ”-g^-] 

♦  [A,(°.  -  b3>  -  »,«»,  -  V]  [Vi(ki>)  ,2 


—ftr*  )  *j"c  v>- 


(19) 


While  the  electric  current  thus  found  is  a  useful  parameter  in  exploring  the  acoustic 
performance  of  a  centrally  supported  disk  radiator,  successful  application  of  four- 
terminal  network  theory  requires  some  simpler  form  than  that  presented  by  Eq.  (19). 

To  obtain  a  more  convenient  expression  for  i ,  we  note,  from  Eq.  (16),  that  to  a  first 
approximation  in  kr ,  the  deflection  w  has  a  parabolic  distribution  with  respect  to  the 
coordinate  r ;  that  is,  W  ~  r2.  If  then  we  assume  a  ueflection  curve  of  second  order  in  r , 

W  =  W0(r2/a2)  . 


we  find 


(second-degree  terms  in  W  only). 


The  charge  for  a  condition  of  parabolic  distribution  becomes 


|Q|  ”  N  »0  +  n*2e*3  E3  , 


(20) 


where  N  =  2ne31  h  =  7ie31 1 . 

The  symbol  N  is  the  ratio  of  electric  charge  accumulating  on  the  electrodes  of  one 
plate  to  the  peak  mechanical  displacement  (at  the  edge  of  the  disk).  To  the  approxima¬ 
tion  of  second-degree  terms  m  r  for  the  variation  of  mechanical  displacement  with  radius 
N  is  seen  to  be  a  real  number.  When  fourth-order  terms  (and  higher)  in  the  deflection 
curve  are  considered,  the  magnitude  of  the  slope  diminishes,  and  with  it  the  transduction 
ratio  N.  To  find  the  transduction  ratio  corresponding  to  higher  order  deflection  curves, 
we  assume  more  complex  displacements,  consonant  with  the  boundary  conditions,  or, 
alternatively,  we  assume  some  known  mechanical  load  on  the  surface  of  the  disk.  Two 
examples  of  the  former  method  and  one  example  of  the  latter  are  given  below. 

Example  1.  Let  the  displacement  curve  be  of  the  form 

This  type  of  variation  with  radius  has  been  used  by  investigp'ors  in  recent  times. 
Southwell  (7)  has  shown  that  such  an  assumed  curve  yields  a  value  of  grave  resonant  fre¬ 
quency  75  percent  of  the  correct  figure,  if  A  -  l  and  b  -  -0.275.  Proceeding  with  these 
assigned  values  of  A  and  B,  we  obtain 
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The  inclusion  of  terms  to  fourth  order  in  the  slope  (dW/dr)  r  -  *  reduces  the  transduc¬ 
tion  ratio  found  for  parabolic  distribution  by  a  factor  of  2.  Correspondingly,  the  presence 
of  an  external  static  pressure  load  introduces  a  similar  reduction  in  n.  Terms  higher 
than  fourth  order  affect  these  results  in  minor  ways  only.  We  conclude  that  the  trans¬ 
duction  ratio  under  actual  operating  conditions  lies  between  In!  =  »r3l  t  and  !n!= (r?/2)e3it. 

We  hare  now  reached  a  point  in  the  analysis  where  further  progress  is  impeded  by  the 
complexity  of  the  approximate  deflection  curve  as  revealed  by  Eqs.  (16)  and  (17).  It  is 
more  convenient,  from  this  point  on,  to  assume  a  deflection  curve  of  3imple  algebraic 
form,  consonant  with  the  boundary  conditions,  such  that  computations  of  acoustic  power, 
etc.,  are  facilitated.  Such  a  choice  would  leave  oqe  factor,  namely  the  peak  displacement 
Vs,  indeterminate.  However,  the  transduction  ratio  N  relates  tc  to  the  charge  Q  and 
therefore  relates  mechanical  force  F  to  applied  voltage  V3.  With  the  magnitude  of  N 
explicitly  known,  the  peek  ;  splacement,  velocity,  etc.,  become  electrical  quantities  whose 
magnitudes  are  then  precisely  determinable.  In  accordance  with  this  procedure,  then,  we 
choose  a  deflection  curve  of  simple  parabolic  form,  namely  W  =  *0  (r/«)2,  and  proceed. 
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Radiated  Acoustic  Power,  Mechanical  Reactive  Power,  Kinetic 
Energy,  Resonant  Frequency,  Acoustic  Pressure  and  Mechanical 
Q  of  a  Centrally  Supported  Disk  Whose  Deflection  Curve  is 
Parabolic  (Second  Order  in  r).  Infinite  Baffle  Present. 

In  the  far  field,  at  a  great  distance  R  from  the  infinite  baffle  in  which  the  disk  with 
parabolic  deflection  is  located,  the  radial  distribution  of  real  acoustic  pressure  p.  has 
been  found  by  McLachlan  (8)  to  be 

pj o*2  2  j2(*)i 

p“ "  ^  L* - (1) 

A- 

*  =  ka  ain  6  k  =  — K  =  wavelength. 

Since  the  liquid-particle  velocity  v  at  great  distances  is  p  ,  the  peak  sound  power 
P,  radiated  into  semi-infinite  space  is 

f  2Jj(z)]2 

P«  =  J  P.  v  dA  =  2*  —  a4 102  J  — 5 - — J  sin  9  d 9 

or 

P  -  <-*>"  Z'W* 

•  '  c«  °  L-*  (2ra+  1)  (2m-  1). . . 

mm  0 

C1-  ,  \2+2m  ft  \4*2«  /,\3«2m 

2 J  +  Mil  f4  +  2m)!  _  4fc)  (3*  2m>!  (22a) 

(2+m)!  (1  +  m)!  (4  ♦»)!  (2+m)!a  (3  +  m)!  (1+m)!  (2  +  m)|J  ' 

When  this  expression  is  expanded,  the  peak  real  power  becomes 

pe  =  2rr  “X2  *4  [is  ~  *  3^6  (k,)4  "  "•]  ‘  (22b) 

At  the  frequency  of  mechanical  resonance  in  the  grave  (umbrella)  mode,  the  magni¬ 
tude  of  ka  is  usually  less  than  1/2.  Limiting  the  infinite  series  of  the  above  equation  to 
the  first  term  only,  we  obtain 


P  "£*.  (peak  power 

•  "  ^  cw  »  ka  «  1/2). 


for  condition 


The  mechanical  reactive  power  Px  can  be  derived  by  a  similar  procedure  applied  to 
the  reactive  pressure  Pi ,  explicitly  written  in  Eq.  (10).  As  in  the  above  case,  we  limit 
the  expression  for  Pj  to  the  first  term  in  ka  and  write 


p*  =  J  Pi  * 

"'A 


=  i  P»c»2naai02  (ka)  J  gjb3*. 


Since  the  value  of  the  integral  is  20/63*,  we  have 


w7V5»,rrrvwr?i«. 
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IS 


P,  =  j  *>L,  W02  (23a) 

«i  =  §  ^.»3  •  (23b) 

The  parameter  Lt  is  the  inertial  mass  which  the  liquid  medium  adds  to  the  disk 
flexing  in  the  grave  mode.  It  coincides  with  the  values  found  by  McLachian  (5). 

The  magnitude  of  the  kinetic  energy  of  vibration  (Tp)  of  the  disk  with  no  acoustic  load 
is  also  found  by  a  simple  integration: 


Tp  =  J  T*rrt  PpW02(^)4  dr 

•V*o 

or 

Tp=|wo2-  M=*»*ppt.  (24) 


The  lumped  mass,  corresponding  to  the  edge  velocity  Wp,  is  therefore  1/3  of  the 
actual  disk  mass.  Adding  to  this  lumped  mass  the  inertial  mass  of  the  water  (Lj),  we 
obtain  the  total  kinetic  mass  Mq,  corresponding  to  wo; 


n  _  120  Pm* 

P  '  63"  p  t  ■ 

P 


(25a) 

(25o) 


The  presence  of  a  semi-infinite  liquid  medium  may  be  thought  to  raise  the-  density  of 
the  plate  from  its  value  pp  to  pp(I  +  /3).  t.  an  ordinary  design  for  a  water  medium  pw«/p  t 
is  close  to  unity,  making  £  about  2/3.  One  result  of  this  added  mass  is  to  lower  the"  P 
natural  resonant  frequency  of  the  disk  in  the  grave  mode  from  its  value  in  a  vacuum 
(approximated  by  air)  by  the  factor  <  1  +  /S )***.  An  expression  for  may  be  obtained  by 
solving  the  secular  equation 

A  =  0  ,  (26) 

where  A  is  the  denominator  of  Eq.  (16)  in  the  absence  of  acoustic  load.  We  obtain  as  the 
first  root  of  this  equation  the  value  ka  =  1.933,  a  magnitude  quite  close  to  ka  =  1.937  found 
by  Southwell  (7).  Upon  solving  for  wQ,  the  vacuum  resonant  frequency  becomes 


Hence  the  resonant  frequency  in  a  liquid  medium  («*)  becomes 


■>*  =  :-««  p  cp  . 


•  _  (_!£_] 

*•  s)J 


(28) 


A  spherical  wave  (such  as  is  radiated  by  the  disk  for  the  condition  ka  «  1/2)  o£ 
source  strength  <?,,  looking  into  a  semi-infinite  medium  bounded  by  a  stiff  baffle,  develops 
a  pressure  p  in  the  far  field  (distance  8)  whose  expression  v4>  J* 


P  = 


-  Q. 


(29a) 
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Now  by  definition,  the  mechanical  Q  is 


Therefore 


“R 

D  «  - —  . 
2ad 


=  2.47(1  *fi)3/2  /fl 


Mechanical  Displacement  Under  Forced  Mechanical  Drive 
and  Receiving  Response 

A  steady-state  (sinusoidal)  pressure  of  magnitude  p0  (independent  of  the  parameter 
ka)  is  allowed  *o  drive  the  submerged,  centrally  supported  disk.  We  stipulate  that 
ka  «  1/2  ana  write  the  equation  of  motion  as 


VI  -  It,  *=$>+«>', 


where 


<p  =  i  <**Z(a,M) 


-  P® 

•  V 


As  for  boundary  conditions,  we  require  that  D3  =  0  (i.e.,  open  circuit).  Assuming 
once  again  that  b/a  «  1,  we  obtain  for  the  first  two  boundary  conditions 

a  +  y  -  <p  +  <p’ 

±0-6  -  0. 

Proceeding  with  the  two  remaining  boundary  conditions,  solving  for  and  6  on 

the  assumption  that  E-  is  zero,  and  substituting  in  the  displacement  equation,  we  obtain 

+!b<)  J  :  A,)] [J0(k ,.) - IoOcj.)] -  [B^A, - Aj) - VB, - M [y0(V>  +1^,.)] 


A  ~  ~ ^  [Bj  (A* (b,  +|b4)J  [Jofk,.)  -  IoOc,.)] 

£VBi-B3)-Bj(A1-AJ))[T(b(k1.)+!K0(k1.)]  -  f^oOcia)-!}  . 


The  complexity  of  this  expression  requires  that  some  simpler  representation  be 
sought  for  the  accumulated  charge  |Q|  due  to  the  pressure  pQ  than  is  contained  in  the 
requirement  that  iQi  =  rgjjhaWdr)^,.  We  assume,  as  before,  a  parabolic  deflection 
curve,  and  calculate  the  open-circuit  voltage  V3  due  to  iQ!  stored  upon  a  capacitor  of 
area 77 »2,  thickness  h,  and  dielectric  constant  e33.  We  have  then 

_V  *31  *2  n 

P.  -  T  .=  <J3  p„  •  I 


where  n  =  2  for  parabolic  displacements.  We  can  use  the  same  formulation  for  the  case 
of  fourth-order  terms  by  setting  n  =  l. 
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It  is  useful,  at  this  point,  to  derive  an  expression  for  the  receiving  response  at  fre¬ 
quencies  so  low  that  Z(a,u)  becomes  negligibly  small,  that  is,  at  quasi-static  load  condi¬ 
tions.  From  the  theory  of  elasticity. 


Po  t3(YB) 

®  AV 


(static  load) , 


where  (Y0  )AV  is  the  average  Young's  modulus  at  zero  current. 
Upon  substitution,  we  obtain 


_V  _  J.  a?  *31n 

p°  4  *  (*<?> 

33  o  'Av 


(static  load). 


(35) 


Formulas  Describing  the  Acoustic  Performance  of  a  Metal-Ceramic 
Bilamellate  at  Mechanical  Resonance,  Radiating«Sound  into  a  Liquid 
Medium.  Parabolic  Deflection  Curve  Assumed  and  Infinite  Baffle  Present. 

At  mechanical  resonance  in  the  medium,  the  magnitude  of  edge  velocity  V0  is  gov¬ 
erned  solely  by  the  radiation  resistance  of  the  medium  to  the  disk  vibrator.  Since  the 
real  acoustic  pressure  for  the  condition  ka  «  1/2  is  independent  of  radius,  the  total 
mechanical  force  restraining  motion  at  mechanical  resonance  is  Fs  =  pn«2.  We  know  that 
edge  velocity  and  pressure  are  related  through  the  impedance  equation,  Eq.  (10).  Hence 

(36) 


When  electrically  driven  by  a  voltage  v3,  the  applied  (mechanical)  voltage  force  Nv3 
balances  the  mechanical  resistance  force  FR  in  the  steady  state.  Recalling  that  N  =  »e31t 
for  the  condition  of  parabolic  distribution,  we  find,  for  edge  velocity. 


or 


4NV 


Pwcw(ka)2  7T  ■ 


2  2  * 
•  TT 


cJE3 

W„  =  l.*l  *31  ~P~ 

CP 


(37) 


This  is  a  key  equation,  for  with  it,  and  with  Eq.  (28)  for  the  resonant  frequency  we 
can  substitute  into  Eq.  (22c)  to  obtain  an  explicit  equation  for  the  acoustic  power 


P. 


5.36  e31 


c_  V, 


Pm 


C 


•2 

P 


(38) 


Similarly,  substituting  Eq.  (37)  into  Eq.  (29b),  we  find  that  the  (real)  acoustic  pressure  in 
the  far  field  (distance  ft)  for  ka  «  1/2  is 


.p!==^^v. 


(39) 


As  for  receiving  response  (v/p,,),  when  an  alternating  pressure  of  constant  magnitude 
p„  is  applied  (at  the  mechanical  resonant  frequency)  to  one  side  of  the  bilamelkte,  and 
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when  simultaneously  the  electric  terminals  are  open-circuited,  an  edge  velocity  is  ini¬ 
tiated  whose  magnitude  is 


4Pq 

P.c.(ka)2  ' 


In  consequence  of  the  electromechanical  coupling  action,  an  alternating  electric  cur¬ 
rent  i=NW0  flows  through  one  plate  of  the  disk  and  charges  the  capacitance  C*  to  a  potential 

_  Q  _  _  *PqN 

"  C*  j  wC*  “  j  P.c.(ka)2^RC*  ' 

Substituting  for  and  N,  we  obtain 


1.582 


*31  c. 


t  e!3p„c 


(40) 


The  open-circuit  voltage  reported  by  this  equation  depends,  for  its  magnitude,  on  the 
presence  of  an  infinite  stiff  baffle.  In  the  absence  of  such  a  baffle,  the  magnitude  of 
response  would  fall  by  a  factor  of  2.  The  response,  in  any  case,  is  that  across  one  of  the 
two  possible  active  plates  of  the  disk. 


Conclusion 

Equations  (6),  (27),  (31),  (37),  (38),  (39),  and  (40)  are  the  most  useful  results  of  this 
part  of  the  analysis.  Taken  to&efer,  they  form  a  relatively  complete  summary  of  the 
acoustic  performance  of  a  centrally  supported  piezoceramic.  One  seemingly  vital  equa¬ 
tion  appears  missing,  namely  an  expression  for  the  motional  admittance.  This  is  not  a 
serious  omission,  sin  e  the  small  magnitude,  of  flexural  electromechanical  coupling 
renders  the  bilamella  almost  a  pure  capacitor,  when  loaded  by  a  liquid  medium.  In 
many  instances  the  m  local  resistance  is  less  than  1/10  of  the  reactance  at  mechanical 
resonance,  and  for  mouc  r.  A  power  absorption  the  current  flow  through  the  radiation 

resistance  is  relatively  sn . 1.  No  great  error  therefore  ensues  in  treating  the  bllamel- 

late  as  a  "pure”  capacitor.  As  for  material  constants,  it  is  best  to  determine  these  by 
actual  tost;  e.g.,  it  is  best  to  determine  cp*  by  measuring  the  resonant  frequency  in  the 
liquid  and  computing  this  quantity  from  Eq.  (28).  Only  in  this  way  can  one  avoid  dubious 
values  for  these  parameters.  Preliminary  estimates  of  the  values  of  these  parameters 
based  upon  the  available  literature  may  be  obtained  from  Table  1,  which  lists  the  piezo¬ 
electric  properties  of  three  popular  piezoceramics. 

An  alternative  procedure  for  determining  the  electrical  series  resistance  R  at  veloc¬ 
ity  resonance  is  to  assume  that  R  %  Xg/Qg,  where  XB  is  the  blocked  series  reactance  at 
resonance  and  Qg  is  the  electrical  Q.  The  factor  Qg  may  be  obtained  from  a  knowledge  of 
and  k,2  by  the  additional  approximation  that  Qg  =  l  -  k  f2/k  f2  At  velocity  resonance 
therefore  the  electrical  Impedance  is  approximately 


z<  <-■«„)  * 
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TABLE  1 

Material  Constants  of  Polarized  Ceramics 


Electromechanical 

Ceramic* 

Property 

(See  List  of  Symbols) 

Ceramic  B 
BaT103  -  CaTi03 

PzT-4t 

PbTi03  -  Pbzr03 

PzT-5t 

PbTi03  -  PbzrO, 

4 

(planar  stiffness) 

127.9*10’  (N/m2) 

91.3  *  10’ 

74.7  *  10’ 

4 

(planar  stiffness) 

38.3  *  10’  (N/mJ) 

27.9  *  10’ 

23.0  *  10’ 

4 

(compliance 

8.62  *10*l2(m2/N) 

12.05  *  10* 12 

12.05  *  lO* 12 

kp 

(planar  coupling 
coefficient) 

0.33 

0.48 

0.54 

c  1 

(strain-electric 
field  modulus) 

-58  *  10* 12  (C/N) 

-97  *  1012 

-140  *  10- 12 

*31 

(planar  piezo  modulus) 

-9.45  (N/mV) 

-10.38 

-12.01 

4  3*3 

(blocked  dielectric 
constant) 

9.64*10*’(C/V,) 

7.89  *  10-’ 

10.24  *  10'’ 

4 

(dielectric  constant  at 
zero  stress) 

10.7  *  io-’(cA.) 

9.7  *  10-’ 

13.22  *  10'9 

(ceramic  density) 

5.4  xlO3  (kg/m3) 

7.6  *  103 

7.6  *  10* 

V 

(Poisson’s  Ratio; 

0.3 

0.3 

0.3 

*The  chemical  content  t*  the  above  listed  ceramics  may  be  obtained  by  writing  tc  the 
author  of  this  paper. 

tThe  two  PbTi03  -  Pbzr03  mixes  are  different  in  composition. 


EDGE-SUPPORTED  BI LAMELLATE  DISKS 

Wa  shall  consider  in  this  part  a  bilamellate  disk  of  two  active  halves,  simply  sup¬ 
ported  at  its  outer  edge.  Care  must  be  taken,  in  building  this  structure,  to  insure  sim¬ 
plicity  of  support,  i.e.,  to  insure  absence  of  clamping  while  at  the  same  time  not  restrain¬ 
ing  radial  displacement  along  the  outer  rim.  In  any  case,  the  boundary  conditions  W  =  0, 

Mr  =  0  must  apply;  that  is,  the  external  casing  roust  restrain  but  not  flex  the  disk. 


Coefficient  of  Electromechanical  Coupling,  Bending  Moments, 

Dynamic  Equation,  and  Solutions 

A  few  of  the  results  of  the  previous  derivations  are  immediately  applicable,  since 
the  mode  of  mechanical  support  does  not  enter  into  the  basic  equations  of  state.  From 
Eqs.  (2a)  and  (2b),  we  obtain  (since  ;  =  2), 

<«*> 

*■»•[»  ~2  t] «»> 
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Similarly,  the  coefficient  of  electromechanical  coupling  in  the  flexural  (grave)  mode  is 


kf 


MIX 


(42) 


In  order  to  write  the  dynamic  equation,  we  need  an  expression  for  p(a,u),  as  required 
by  Eq.  (9).  Now  the  radial  distribution  of  pressure  p(r.t)  has  a  form  identical  with  that 
expressed  by  Eq.  (10),  with  the  exception  that  the  g's  (hypergeometric  functions)  contain 
additional  terms.  Restricting  ourselves  once  again  to  the  condition  ka  «  1/2,  we  employ 
only  faciors  and  g2  in  the  infinite  series  of  Eq.  (10).  On  the  assumption  that  the 
deflection  curve  has  the  form  w  =  wc[i  -  (r/a)2] ,  we  consult  Me  Lachlan  (S)  and  find  that 
g2  =  1/2,  and  is 


>•  bJ)  +  b*  [F  (" 


I  I  X 
2’  2’  * 


Performing  once'  again  the  necessary  integrations,  we  find  that  the  reaction  pressure, 
averaged  over  the  surface  area  of  the  bllamellaie  disk,  is 


p(a,u)  =  j  ulc  Z *(a,w)  (43a) 

Z'(a.u)  =  pmcw  +  j  ^  ka  j  .  (43b) 

It  may  be  concluded  that  a  change  in  the  mode  of  support  has  altered  the  pressure 
distribution  in  a  minor  way  only,  assuming,  as  has  been  done,  that  the  deflection  curve  is 
parabolic  in  both  conditions  of  support.  Equations  (43a)  and  (43b)  lead  directly,  in  con¬ 
junction  with  Eq.  (13),  to  the  equation  of  motion  (steady  state): 


4  .  j  U»,  Z'(*,u) 

y  •  -  V4 V  =  - _ _ _ _ _ 

1  D* 


(44) 


Although  Eq.  (14),  with  replaced  by  Wc,  is  a  general  solution  of  the  differential 
equation,  the  boundary  conditions  noted  in  the  introduction  to  this  section  of  the  analysis 
require  that  6  and  b  be  zero.  Our  general  solution  to  Eq.  (44)  is  therefore 


W(r,u0 


*  Jo<klr)  +  Tl0(k|O  ~ 


j  uZ'(*,w) 

D*  k4 


{45) 


The  two  conditions  that  determine  a,  y  are 


aJo<ki*)  ♦  r!o(ki«)  = 


D*  k  4 
l 


a*l  +  >*3  = 


*31  S3 

2  0  kJh 

11  l  n 


Upon  solving  these  simultaneously  and  substituting  the  results  in  Eq.  (45),  we  obtain 


3  *31  ®.l  J"lofkI*5  Jo<ki«->  *  JoOti*)  I0(kif) 

2  i°~k)h  L  a7 


.  j  u'  Z'(m,u>)  J"43  Jofkir)  ~  Io(kjr)"J  j  uZ’(n,u) 

D*k4  *L  A*  J-  D-k7“fc  . 
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where  A'  =  A3  J0(kj«)  -  Ai  I0(kja) . 

From  this,  upon  setting  r  =  o,  we  obtain 


”31  t3 


11 


_ fI0(M)  -  J0(M)3 _ 

t  ju>  A'Z'(a,o>)  j£i)2'(a,o<)(A3 -Aj) 

A  +  D*k4  D*  k 4 


(46) 


(47) 


The  displacement  V,  according  to  Eq.  (46),  is  parabolic  to  a  first  approximation  in 
the  parameter  kr.  We  have  seen  (Eq.  20)  that  for  such  a  distribution  of  deflection,  the 
transduction  ratio  N  for  one  plate  of  the  disk  is  ^e31 1.  Under  surface  loads,  however, 
the  transduction  ratio  diminishes,  as  may  be  seen  from  the  following  development.  Let 
the  simply  supported  disk  be  subject  to  a  static  surface  pressure  pq.  The  charge  accu¬ 
mulating  on  one  pair  of  plate  electrodes  is 


Q  = 


a31 


r 


♦  *t 
2 


2n  r  dr  , 


where 

d31  is  the  electric  displacement-stress  "effective"  piezo  modulus 
sr.st  are  the  radial  and  tangential  stresses  induced  in  the  disk  by  the  load  pQ. 
From  the  theory  of  thin  plates,  it  may  be  found  that 


and 


+  s. 


O  C 


Upon  performing  the  required  integration,  we  obtain 

<?  =  J"Vd31  *= 


w*31*  *c  • 


that  is,  the  transduction  ratio  N  for  the  (static)  loaded  state  is  3/4  of  the  value  predicted 
by  a  parabolic  deflection  of  the  disk.  When  two  plates  are  considered,  the  ratio,  of 
course,  will  be  twice  the  value  of  one  plate.  The  advantage  of  using  two  active  plates  is 
thus  clearly  evident,  since  the  coupling,  power,  pressure,  etc.,  will  all  be  improved. 


Acoustic  Power,  Mechanical  Reactive  Power,  Kinetic  Energy, 

Mechanical  Resonant  Frequency  and  Mechanical  Q  of  a  Simply 
Supported  Disk  Whose  Deflection  Curve  is  Parabolic  (Second 
Order  in  r).  Infinite  Stiff  Baffle  Present. 

An  expression  for  the  acoustic  pressure  in  the  far  field  'distance  R)  due  to  a 
baffled  disk  having  a  deflection  curve  of  the  form  W  =  *cfl  -  (r/*)7]-  is  given  by 
Me  Lachlan  (5): 
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2p„«JWe  J2(kja  sin  9) 

P.  =  — T—  (kia  sin  e)2  •  (48) 

A  suitable  integration  over  the  spherical  area  of  radius  R  of  the  product  of  the  acous¬ 
tic  pressure  and  the  particle  velocity  yields  the  peak  acoustic  power  radiated.  symbols, 

p.  =  P.vdA  =  I  P2rfA 


w  -2  4  f1  h  (k!a  »in  g>  . 


o  (kj»  sin  9)4 


Upon  evaluating  the  definite  integral,  we  obtain 


sin  6  A6. 


/i\2(2*«)  - 

4-2  y1  (2)  (k>)  m2aml 

•  n  cw  *  c  Z_j  m!(4  +  m)!  (2  +  n)!J  (2m  +  1)(2»- 1)  . . .  * 

■C  0 

or 

We  note  that  if  ka  «  1/2,  the  expression  for  peak  acoustic  power  in  the  medium  due 
to  a  flexible  simply  supported  disk  is  the  same  as  for  a  centrally  supported  disk,  both 
having  parabolic  velocity  distribution  and  both  having  equal  maximum  displacements. 

To  terms  of  first  order  in  ka,  the  reactive  pressure  on  the  surface  of  the  disk  is 

Pi  =  i«Vw*c(ka)gj. 

The  magnitude  of  reactive  power  in  the  medium  (Px)  is  therefore 

f 

P,  »  jw277  a3pwWci  |  gj(l-b2)  bdb 
0 

or 

px=il  -P„a3*c2  -  (50) 

The  mechanical  reactive  power  of  the  medium  whose  expression  has  just  been 
derived  is  due  entirely  to  the  inertial  effect  of  the  medium.  If  the  motion  of  each  ele¬ 
mentary  volume  of  adjacent  liquid  is  referred  to  the  peak  velocity  Wc,  we  see  that  the 
inertial  mass  added  by  the  medium  is  p„«x3  (248/315).  Since  the  kinetic  energy  of  the 
vibrating  plate  Vs 

# 

Tp  =  |J  PwcwWc2(ka)27T«2(l-b2)bdb 
0 

r  5  »  «2  V  *<?  - 

we  conclude  that  the  total  effective  mass  of  plate  and  water  (u ')  referred  to  the  peak 
velocity  ic  is 


A 

J 
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M;  =  iiraJppt(W) 

^(Sf)^'0-757^-  (5« 

By  thus  increasing  the  effective  mass  of  the  plate,  the  liquid  medium  acts  to  reduce 
the  natural  frequency  of  vibration  in  air  by  the  factor  (1  +  /5')*H  Ah  expression  for 
a/  is  obtainable  from  Eq.  (46)  by  setting  Z(a,u>)  equal  to  zero  and  solving  the  secular 
equation  A'  =  o  for  the  lowest  root,  excluding  ka  =  o.  We  obtain,  for  this  lowest  root,  a 
valve  ka  %  2.252,  from  which  we  find  that  the  grave  resonant  frequency  of  a  liquid-loaded, 
simply  supported  plate  to  be 

uL  =  1.468  -4  c'* 


\ppa  +  (8*T 


Since  for  equal  maximum  displacements  (WG  =  wc)  the  radiated  acoustic  power  is  the 
same  whether  the  disk  is  simply  supported  or  centrally  supported,  and  since  for  parabolic 
displacements  the  kinetic  energies  of  the  unloaded  plate  are  also  identical,  we  see  that  the 
temporal  damping  factor  has  the  same  form  as  in  Eq.  (30),  with  the  exception  that  /S  is 
replaced  by  /9'  and  <^R  is  replaced  by  The  mechanical  Q,  therefore,  becomes 


3/2  pcm 


where 


Cp  = 


Formulas  Describing  the  Acoustic  Performance  of  a  Simply  Supported 
Ceramic  Bilamellate,  at  Mechanical  Resonance,  Radiating  Sound  into  a 
Liquid  Medium.  Parabolic  Deflection  Curve  Assumed,  and  Infinite 
Baffle  Present. 

In  this  section  the  same  reasoning  and  the  same  steps  in  derivation  are  applied  as 
was  used  in  obtaining  Eqs.  (36)  to  (40),  with  the  additional  consideration,  however,  that 
the  transduction  ratio  N  is,  for  the  entire  disk,  2X  the  value  previously  used.  A  first 
result,  employing  is  the  center  velocity,  wc.  We  obtain,  for  this  key  parameter. 


W,  =  1. 85 


531c»E3 


Similarly,  using  Eqs.  (19)  and  (52)  and  this  value  for  wc,  we  derive  the  expression 
for  the  acoustic  power  radiated 


2  C«V3 

r.  =  »-«  *3i  - — ttj 

PwCr. 


As  for  the  acoustic  pressure,  it  is  noted  that  when  z  «  1/2  the  limit  of  J2(z)/z*  is 
1/8.  Hence  the  resonaht  acoustic  pressure  in  the  far  field  (distance  R)  on  the  acoustic 
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axes  has  the  same  form  as  in  Eq.  (29b),  with  the  exception  that  Wc  replaces  wo  and  <4 
replaces  ^ .  Substituting  as  before  for  -4  and  wc,  we  obtain 

lp|  =  1.36  e31  .  (56) 

c'  H 

A  procedure  similar  to  that  which  led  to  Eq.  (40)  yields  the  resonant  open-circuit 
receiving  response 


0.632 


e 


31 

Pm 


(57) 


We  observe  here  the  same  precaution  in  the  application  of  Eq.  (57)  as  was  noted  in 
the  application  of  Eq.  (40),  namely  that  this  is  the  response  from  a  single  plate  (of  two 
possible  active  plates)  and  that  an  infinite  stiff  baffle  bounds  the  half  space  from  which 
the  incoming  signal  takes  its  origin.  In  the  absence  of  the  baffle,  the  response  will  fall 
to  a  value  of  one-half  that  noted  above. 


Conclusion 

Equations  (49),  (52),  (53),  (54),  (55),  (56),  and  (57)  constitute  in  their  entirety  a  sum¬ 
mary  of  the  acoustic  performance  of  a  simply  supported  flexural  bilamellate  disk  radiat¬ 
ing  sound  into  an  infinite  half  space.  All  the  limitations,  precautions,  etc.,  noted  in  the 
discussion  at  the  conclusion  of  the  earlier  derivations  are  applicable  to  these  last  results. 
In  particular,  it  is  of  importance  to  repeat  the  stipulation  mentioned  previously  that  all 
material  constants  occurring  in  the  formulas  be  determined  by  test  upon  an  actual  disk. 
Substitutions  from  generalized  data  available  in  the  literature  may,  or  may  not,  lead  to 
dependable  results.  And  a  final  point:  all  power,  pressure,  and  voltage  response  for¬ 
mulas  are  based  on  assumed  100-percent  energy  conversion,  no  losses  occurring  on  the 
way.  In  actual  practice,  overall  conversion  efficiencies  vary  from  40  to  70  percent  at 
low  electric  drive  (0.01  volt  rms  per  mil  of  thickness)  to  15  to  20  percent  at  high,  elec¬ 
tric  drive  (10  volts  rms  per  mil  of  thickness).  In  practice,  too,  the  material  constants 
e31  and  eft  are  lossy  (i.e.,  are  complex  quantities)  and  are  frequency  sensitive.  From 
these  remarks,  the  approximate  nature  of  the  derived  equations  may  be  surmised. 
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LIST  OF  SYMBOLS 

Symbol  Dei’lnition 


A 

a 

C* 

cn,  c12,  e:c. 

c. 

D 


f 

r1,  g3,  etc. 
li 

J„(z).  In(z). 

v*>.  VO 

k 

kl 

kf 

Mr 

M, 

N 

P. 

P« 

Pi 

0 

<in 


area  in  square  meter 

radius  of  plate  (meter) 

capacitance  (farads)  at  constant  strain 

stiffness  moduli  (N/m2) 

velocity  of  sound  in  water  (m/sec) 

3tiffness  constant,  (N-m) 
electric  displacement  (coul/m2) 
piezo  modulus 

electric  field  (volt/m1) 
frequency  (cps) 
hypergeometric  functions 
hall  thickness  of  plate  (m) 

Bessel  Functions  of  order  n  and  argument  z 

wave  number,  (m*1) 

3  Pp  a.2/h2  cvDi 

coefficient  of  electromechanical  coupling  in  flexural  mode 

radial  bending  moment  per  unit  of  length  (Nm/ra ) 

tangential  bending  moment  per  unit  of  length  (Nm/m) 

transduction  ratio 

acoustic  power  (watts) 

real  pressure 

inertial  reactive  pressure 

charge  (coulombs) 

mechanical  Q 
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Symbol 


r.  B,  i 


a.  fS,  y ,  b,  r 
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LIST  OF  SYMBOLS  (Continued) 

Definition 

internal  resisting  shear  per  unit  of  circumferential  length 

distance  in  sound  field  (m) 
polar  coordinates 
radial  strain  (m  n) 
tangential  strain  (r:  h) 
radial  stress  (N  -2) 
tangential  stress  (N  r3) 
applied  voltage  (volts) 
velocity  in  meters  per  second 

plate  deflection  in  2  direction  as  a  function  of  time  (m) 
plate  deflection  in  2  direction  as  a  function  of  frequency  (m) 
subscripts  for  water,  plate 

Young's  modulus  (N.  m2) 
acous*;c  impedance 
k  n 

constan1  - 

dielect -ic  ( yistant  at  constant  (i.e.,  zero)  strain  (farad/m) 

1  or  2,  depending  upon  backing  plate  of  bilamellate 
density  <V_  r ;) 

Poisson's  ratio 

angular  frequency 

special  denominator  in  Eq.  (16) 


NOTE:  The  MKS  system  is  used  throughout  this  paper. 
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